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ABSTRACT 



I present a matched- wave asymptotic analysis of the driving of solar oscillations by 
■ a general localised source. The analysis provides a simple mathematical description 

. of the asymmetric peaks in the power spectrum in terms of the relative locations of 

eigenmodes and troughs in the spectral response. It is suggested that the difference 
\ in measured phase function between the modes and the troughs in the spectrum will 

provide a key diagnostic of the source of the oscillations. I also suggest a form for the 
| asymmetric line profiles to be used in the fitting of solar power spectra. 

Finally I present a comparison between the numerical and asymptotic descriptions 
of the oscillations. The numerical results bear out the qualitative features suggested 
by the asymptotic analysis but suggest that numerical calculations of the locations of 
the troughs will be necessary for a quantitative comparison with the observations. 



Oh 

2 ' Subject headings: Sun: oscillations 



1. Introduction 

Helioseismology is primarily concerned with the peaks in the solar oscillation spectrum which 
are identified with acoustic eigenmodes of the Sun. However the whole spectrum, including the 
variation between the mode peaks, contains useful information about the Sun's structure and 
dynamics. Moreover, the asymmetries in the p-mode lines which were first identified by Duvall 
et al. (1993) can bias the fitting of profiles to the peaks and hence introduce errors in the 
determination of the frequencies. The most recent space-based spectra from SOHO-MDI (Nigam 
et al. 1998) have confirmed the existence of a strong asymmetry. Thus the study of p-mode line 
asymmetries is prompted both by the need to improve on traditional data analysis techniques and 
by the hope of obtaining new scientific information from them. 

The occurrence of highly asymmetric line profiles is a common feature of models in which 
the oscillations are driven by a localised source (Gabriel 1992, 1993, 1995). Abrams and Kumar 
(1996) explored the nature of the asymmetry in a simple model using analytical techniques and 
then applied these ideas to the interpretation of numerical calculations of the line asymmetry in a 
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real solar model, including the variation of the asymmetry with frequency and mode degree and 
estimates of the error introduced into the determination of the frequencies by the neglect of the 
line asymmetry. 

Rast & Bogdan (1997) emphasised the sensitivity of the asymmetry to the depth and type 
of the source. In particular they introduced the approach that the observed asymmetry can 
best be approached through understanding the relative locations of the peaks and troughs in 
the spectrum, where the location of the troughs can be determined (in some circumstances) by 
solving an eigenproblem for a solar model truncated from above at the source depth. Although 
the calculations presented here will result in an estimate for the whole spectrum, I will follow the 
example of Rast & Bogdan in concentrating particularly on the relative locations of the peaks and 
troughs as a key diagnostic indicator of the source of excitation. 

Roxburgh & Vorontsov (1995) used a form of asymptotic analysis to model the power 
spectrum. The analysis presented here differs from theirs in that it uses a matched wave technique 
in which the matching across the turning point(s) is treated explicitly. Moreover, I consider a 
general form of the source term and obtain an explicit expression for the asymptotic line profile. 
It is found that the power spectrum can be represented in a simple form even for a general 
source, and this should guide the analysis and interpretation of the peaks and troughs in the solar 
p-mode spectrum. 

I begin in §^ with a very simple illustrative model which nevertheless captures some important 
aspects of the dynamics. In §|| I show how one may carry out a matched-wave asymptotic 
analysis for a stratified layer with a source containing both monopole and dipole terms. I am 
predominantly interested in the lowest frequency modes, for which the asymmetry is most easily 
observed. However, asymptotic results for intermediate and high-frequency oscillations are also 
presented. In §|4] I evaluate the resulting asymptotic expressions for a simple model and compare 
them with a direct numerical calculation. 



Figure [T] shows a simple model of a driven oscillator consisting of a uniform string characterised 
by a soundspeed c, clamped at its ends at z = 0, L and driven by a localised source at z = zq <C L. 
The displacement, of the string satisfies the equation 



where the wavenumber k = uj/c and the amplitudes of the monopole and dipole source terms, ao 
and a%, are, in general, arbitrary complex functions of the frequency lo. 

Although such a model is almost trivial to analyse, its solution will be found to illustrate 
many of the properties of the more general problem with which we are primarily concerned. We 
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An Illustrative Model 




(1) 
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can obtain jump conditions at z = zq by integrating the equation of motion from z = zq — e to 
z = zq + e and by multiplying the equation by z and then performing the same integral: 



dl 



d^ 

cL7 



= a (2) 

z -e 



zo+e 

and 

*Zo+e ~ *«d- 6 = ax. (3) 

Evidently a dipole source is inconsistent with a physically realisable string, which must be 
everywhere continuous, but it is a useful idealisation; equivalent to a pair of oppositely-signed 
monopole sources whose separation is much less than any other length-scale in the problem. The 
solution to equation @j is then 

9= (A B mkz, z<z , 
\ B sin k(z — L) z > Zq, 

where 

A = - — \—kai cos k(zo — L) + an sin k(zo — L)\ (5) 

k sin kL 

and 

B = — [— kax cos kzn + an sin kzn \ . (6) 

ksmkL 

It is immediately evident that the amplitudes both diverge for kL = nir, i.e. at an eigenmode of 
the whole cavity. 

Now in the Sun, the driving occurs close to the surface but below the region where the 
observations are made so, by analogy, we define the measured power as \A\ 2 . However the total 
energy in the string is proportional to zn|^4| 2 + (L — -Zn)|-B| 2 . For convenience we define a quantity 

ao + ax ,_s 

/i = . (7) 

a - ax 

Figure || shows plots of the power and energy spectra for several choices of \i and zq . The following 
properties of the solutions are clearly evident: 



(i) The degree of asymmetry in the power spectrum increases as the source is moved closer to the 

upper boundary. 

(ii) The peaks in the energy spectrum are much more symmetric than those in the power 

spectrum. Thus, in a sense, the zeros in the power spectrum are illusory, or at least 
misleading, since there are no zeros in the energy spectrum excited by the source. 

(iii) For the lowest frequency modes, and for a monopolar source, the troughs lie above the 
frequency of the corresponding modes (a negative asymmetry in the language of Abrams & 
Kumar). 
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(iv) For a monopole source, the asymmetry decreases as the frequency is increased and eventually 
changes sign. 

(v) A purely dipolar source results in peaks with the opposite sense of asymmetry to that given 

by a monopole source. 

(vi) For a dipole source, the asymmetry increases as the frequency is increased and eventually 
reverses sign. 

As we shall show below, all these properties are also present in the power spectrum for a 
stratified medium with the exception of property (vi). Instead, for the stratified case, a purely 
dipolar source located in the evanescent region gives an asymmetry which decreases in magnitude 
with increasing frequency. 



3. Asymptotic Analysis for a Stratified Atmosphere 

We restrict ourselves to consideration of the oscillations of a plane-parallel stratified layer and 
ignore the effects of curvature. Oscillations of such a region with frequency to and horizontal wave 
number k satisfy the equation 

+ K 2 * = A^z) (8) 

where the independent variable z increases in a downward direction and the dependent variable 
^ = p x / 2 c 2 V • u, where p and c are the equilibrium density and sound-speed respectively and u is 
the velocity field of the perturbation. The form of the source term, A^^z), is discussed in detail 
by Gabriel (1993). The local vertical wavenumber, K, satisfies 



2 
C Z 



where 

„2 



^ _ J2 _ c 2 k 2 1 _ 11_ 



(9) 



N 2 = g (i^P_J_d}np\ (u) 



and 

N 2 = a('- 

dz Ti dz 

and where p and g are, respectively, the equilibrium pressure and the acceleration due to gravity 
and T\ = {d\np/d\n.p) S: the derivative being at constant entropy. 



3.1. Driving in the Evanescent Region 



For a solar-like stratification, and for frequencies less than the photospheric cutoff at about 
5.3mHz, there is a cavity in which K 2 > 0, bounded by turning points at which K 2 becomes 
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negative. The lower turning point is dictated by the rapid increase of c with depth, while the 
upper turning point is dictated by the rapid increase of oj c near the solar photosphere. Modes 
close to the cutoff have an upper turning point very close to the surface while lower frequency 
modes are reflected somewhat more deeply. 

The precise form of the source terms in equation (||) cannot be deduced without a complete 
theory of solar convection. However, it is clear that the source must be strongly peaked in the 
turbulent super-adiabatic boundary layer at the top of the convection zone, a result also evident 
from the observationally-determined excitation rate (Goldreich et al. 1994). In this section we 
will restrict our analysis to the case of low-frequency oscillations for which this source lies in the 
evanescent region above the upper turning point. 

For simplicity we will consider the case of a source localised at a single height: 



where the factors |.ff(zo)| 2 have been introduced purely for convenience (since ao and a\ may, 
in general, be functions of frequency and wavenumber). The solution for a more general source 



ai as monopole and dipole terms respectively. This use of the terms "monopole" and "dipole" is 
quite distinct from that employed by Goldreich & Kumar (1990) who define monopole terms as 
those sources involving a volume change, dipole terms as those arising from a momentum source, 
and quadropole terms as those arising from internal stresses. This latter definition is clearly the 
better physical description of the actual excitation mechanism. Moreover, as Goldreich & Kumar 
point out, a definition of the multipole expansion solely in terms of the number of derivatives 
appearing in the source terms is ambiguous unless one also specifies ones preferred choice of 
independent wave variable. However, given a particular choice of wave variable (and the variable 
^ is a natural choice) the form of the power spectrum is indeed determined by the number of 
derivatives (or rather the number of sign changes) appearing in the source terms. The precise 
relationship between the two descriptions of the source terms depends on the detailed properties 
of the turbulent convection itself. In the future we hope to expand on the numerical work of Stein 
& Nordlund (1991) to obtain a fully consistent picture of the excitation process. However, for 
the present purposes, the terms "monopole" and "dipole" will be used to refer to terms involving 
ao and a\ as above, with no particular assumption being made as to the physical nature of the 
driving mechanism. 

The geometry of the problem is indicated in figure (|3|). The propagating region (Region II) 
is bounded by the lower and upper turning points z% and Zt respectively, and the source at zq 
divides the upper evanescent region in two. We therefore proceed by writing down approximate 
solutions to equation (||) valid in each of the four regions and using matching conditions at the 
boundaries between the regions to fix the arbitrary constants appearing in the solution. 



A,,k(z) = a \K(z )\U(z - zq) + ai\K(z )\ U'(z-z ) 




term will be given in 



j |3. 1.2 . In the following analysis I shall refer to the terms involving ao and 



I shall use the (lowest order) WKB solutions, which are valid provided the scale of variation 
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of K 2 is long compared to the vertical wavelength, i.e. 



d\K\ 



dz 



< \K 2 



(13) 



This condition will certainly break down close to the turning points where K = 0. The lowest-order 
WKB solutions for the geometry of figure |3| are 



I : * ~ 

II : f ~ 

III : * ~ 

IV : * ~ 



— — r^4exp 
\K 2 \1 



r zt \K 2 {z>)\i&z> 



-±r {B exp \i J* K(z')dz'] + C exp [-i J£ A" (s')de'] } 



\K 2 \? 



Dexp 



Ji K 2 (z')l*dz' 



'20 

1 

\K 2 \? 



+ -Eexp 



/' \K 2 (z')\*dz 



Fexp 



2 » Zt 

z t » z > Z t 

Zj > z > z 

z < z (14) 



where ^4 — F are (as yet undetermined) constants. A requirement for the validity of this analysis 
is clearly that Zq <C Zf, which is to say that the source must lie in a region in which the waves are 
highly evanescent - the low- frequency limit. 



The jump conditions at z = zq are determined exactly as in §H 



d^ 

d7 



dm 

~dz~ 



a |iT(zo)| 2 



and 



giving 



and 



z -e 



ai\K(za) 
-D-E + F = a 1 
D - E + F = -ao 



(15) 

(16) 
(17) 
(18) 



when derivatives of K 2 are neglected. 



The determination of jump conditions at the turning points is essentially standard analysis. 
In the region of the upper turning point, K 2 is assumed to vary linearly so that 



K 2 = P(z - Z t ) 

and the solution can be written 

* « aiM(p*(Z t -z))+ a 2 B[(p l 3(Z t - z)) 



(19) 



(20) 



where Ai and Bi are Airy functions and a\ and a 2 are constants whose linear relationship to 
B,C,D and E is to be determined. We can substitute the expansion ( |l9|) into the WKB solution 
(14) to obtain 

« /H(z - Z t )-\ (fl e §tf* + Ce-^*-^\ (21) 
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in the region z > Zt and 



(22) 



in the region z < Z±. Here we have introduced 



A 



K 2 \2dz. 



(23) 



The Airy function expansion (eq. |20| ) can be expanded asymptotically for large negative or 
positive arguments to give 



$ Ft 7T 2/3 12 (z — Zt) ] <j <\ I Mil 

and 



2 1 , ,3 7T 



+ Q2 COS 



2 1 , ,3 7T 



z > Z t 
(24) 

(25) 



It is now straightforward to match the coefficients of exponentials in the overlap region 
between equations (pi]) and ( |24| ) and between ( |22"|) and (^5|) to obtain four equations from which 
one may eliminate a\ and 02 to obtain 



BeT~ + Ce~ = 2iEe l 



and 



Similar analysis at the lower turning point gives 
where 



Kdz. 



(26) 
(27) 

(28) 
(29) 



The five equations (17, 18, 26, £7] and 28) determine the five amplitude coefficients B-F. (A 
second condition at the lower turning point can be used to determine A if it is required.) 

In the solar case, the oscillations are measured above both the upper turning point and the 
source so the measured amplitude spectrum will be the quantity F (multiplied by some slowly 
varying function of frequency to convert the spectrum of Vl/ to a doppler velocity spectrum). It is 
straightforward to solve the equations to obtain 



1 



1 



F = -(-ao - 01) + 7 (-00 + ai)e 



-2A 



tan ( 



(30) 



which is the low-frequency asymptotic amplitude spectrum. 
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3.1.1. Properties of the Power Spectrum 

The spectrum denned by equation (|3^) diverges at frequencies corresponding to eigenmodes 
of the cavity, i. e. 

' Zt Kdz= (n- 7T (31) 



independent of the form and location of the driving. However the minima of \F\ 2 depend on the 
values of Zo, ao and a\. Thus, as we would expect, the spectral response away from the eigenmodes 
contains information about both the cavity and the source. If we consider cases where the ratio 
(ai — ao)/(ai + ao) 1S rea l then the minima will, in fact, be zeros (nulls). We are particularly 
interested in the low- frequency /large- A case, for which it is evident that the zeros must lie close 
to the eigenmodes, i.e. for which the peaks in the power spectrum must be most asymmetric. In 
this case it is easy to show that the zeros are given by 

rKdz*(n- 1 -)n + 2e-^^±^. (32) 
Jz t V 2/ a -ai 



Now the integral in equation ( |29| ) is known (Christensen-Dalsgaard & Perez Hernandez 1992) to 
have the property 



Zt 



Kdz^u;F D +I(u) (33) 

where the function Fu(uj/k) is the Duvall function and I(uj) is a term whose properties depend 
predominantly on the structure of the surface layers. Furthermore, it is clear that A is almost 
independent of k since it involves the evaluation of K only in the evanescent region where u S> ck. 
Thus both the dispersion relation [eq. [HJ and the equation for the nulls [eq. [3^] can be written in 
the form 

Mih*— - (34> 

where the Duvall function is the same in each case but the phase function a(uj) for the troughs 
differs from that from the modes. In the limit of large A we have 

at ~ a m = -e (35) 

7T ao — CL\ 

where a t and a m are the phase functions specifying the positions of the troughs and modes 
respectively. The sign of the difference depends on the type of the source. For a pure monopole 
source, at — a m is positive, corresponding to troughs in the power spectrum lying above the 
corresponding peaks, with the situation being reversed for a pure dipole source. In §|| I will show 
results for a simple model which illustrate these properties of the power spectrum. 

The special case ao = a\ corresponds to a solution with A = B = C = D = E = 0. In 
this case there is no wave energy below the source and so the spectrum contains no signal of the 
presence of an acoustic cavity. Thus the spectrum of F is flat. 



- 9 - 



3.1.2. The General Solution 



If we consider the solution to the local-driving problem [eq. [T^] with a\ set to zero then we 
immediately recover the Green's Function for the general case of equation (||). Thus, if the source 
term is located entirely within the evanescent region one can simply write down the solution in 
the region above the source: 



exp 



JJ I XV) ids' 



A\K 2 {z)\± 



J A^zo) {2e A ^ +e- A ^Han^} \K(z )\-^ dz 



(36) 



where A and (f> are defined by equations (23) and ( |29D respectively. Note that 4> is independent of 
Zq. Thus multiplying equation (|36| ) with its complex conjugate one obtains the power spectrum: 



W ■ l\IU:)\ 2 \^Z\ 2 h + 9?(|£Jtan0+i 



s: 



tan 



(37) 



where I have introduced 



and 



A 



LU,k( Z 0) e 



±A( 20 ) \ K[ 



i 

~ 2 dzn 



exp 



Hjz) 



fl \K 2 (z>)\idz> 



(38) 



(39) 



4|K 2 (z)|3 

Treating equation (^) as a quadratic in tan^>, it is clear that real zeros (true nulls of the 
spectrum) can exist only if S~ / is real. Thus if we were to replace our delta- function source 
by a finite-width, single-phase source, the troughs in the spectrum would still be true zeros. 
Conversely, the existence of sharp troughs in the measured spectrum cannot be taken as an 
indication that the source itself is necessarily highly localised. 



3.1.3. Proposed Line Profile for Low- Frequency Modes 

A key question in helioseismology is the correct way to determine the true solar oscillation 
eigenfrequencies from the measured spectrum of asymmetric, multiply-excited peaks superimposed 
on a noisy background. Rast & Bogdan (1997) suggested that a correction to the peak location 
might be obtained by using the location of the nearest trough while Rhodes et al. (1997) 
show a comparison of an asymmetric p-mode line profile to an asymmetric Lorentzian function. 
The asymptotic spectrum (|^) suggests a possible physically-motivated profile for use in fitting 
asymmetric peaks in the data. For a given peak one takes 

< f >= I( y - Vo ) + l + ia 1 (40) 
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where S is a typical inter-peak spacing and the damping constant a\ has been introduced to ensure 
that the amplitude is everywhere finite. Then the peak profile will be given by 



\F\ 2 = a 2 + «3 



- — •■■ - 1 + (a 4 + ia 5 ) tan ( %{y — uq) + ^ + iai 

o 2 



2 



(41) 



where a\-a^ and (5 are real parameters. The additive noise background has been expressed as a 
constant a 2 , although in practice a power law might be preferred. This can be written 

|-F| 2 = a 2 + a 3 x 

{(o!4 cos x — sinx) cosh a\ — a$ sin x sinh ai} 2 + {(04 sinx + cosx) sinh a\ + 05 cos x coshai} 2 

sin 2 x cosh 2 a\ + cos 2 x sinh 2 a\ 

(42) 

where 

x = ^{u-u ). (43) 



One could fit the peak directly to (42). Alternatively one notes that close to the peak we have 
tan <f> pa —1/(0 — vr/2) so 

- -c n 2 - u :! , 

x z + a\ 

Equation p^ ) produces only a single asymmetric peak while equation ( f42| ) includes the influence 
of nearby peaks directly due to the periodicity of the tan function. Thus it is clear by inspection 



\F\^a 2 + a 3 {x - Cu)2 2 +ia r a5)2 (44) 

t<5 -J- 



that equation (44) cannot be used to fit the parameter 5 independently of the other parameters. 



An important limit of equation (|44| ) is the case 0:5 = which yields 



\F\^a 2 + a 3 + 4^(l-^) (45) 
x l + af \ 04/ 



from which one can see that allowing 04 to tend to infinity while keeping a\a 3 constant yields 
a symmetric Lorentzian peak. Thus the asymmetric forms include symmetric peaks as a special 
case. 



The actual use of either (42) or ( (44| ) to fit solar data will be the subject of a future 
investigation. However, in order to illustrate that such a procedure is feasible, figure |3| shows some 
preliminary results of a fit of profile (|4^) to one month of low degree GONG data. Even from a 
single month of data, it is clear that the peaks show a negative asymmetry (more power in the 
low-frequency wing of the peak than in the high-frequency wing) and that the asymmetric profile 
produces a substantially better fit to the data than does a symmetric Lorentzian. A full discussion 
of the fitting procedure, error analysis, and details of the numerical results will be presented in 
a future publication, and will be based on analysis of much longer observational runs from both 
GONG and the SOI-MDI instrument on SOHO. 
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3.2. The High-Prequency Case 



In the Sun, oscillations are seen well above the photospheric cutoff frequency. These have 
been interpreted ( Kumar et al. 1990, Kumar & Lu 1991) as pseudomodes (interference fringes) 
caused by the interference of a wave propagating up from the source with a wave propagating 
downwards from the source and reflected back upwards at its lower turning point. The analysis 
developed in § [3,1| is easily extended to this case. We now have only three regions to consider: a 
lower evanescent region, a central propagating region below the source and an upper propagating 
region above the source. The WKB approximations to the solutions in these three regions are 



I : 

II : 
III : 



— r ^4exp 

K 2 \Tt 



r zt \K\ Z ')\-^z> , 

V {B exp [t fz t K(z')dz'} + C exp [-« J| ( K(z')dz'] } , 

is: K(z')dz' . 



-^-D exp 

K7. 



z > z t 

z t » z> z 

Z < Z 



(46) 



Using the same matching condition at z% and jump conditions at zq leads to the amplitude 
spectrum in the upper region: 



where 



Kdz. 



(47) 



(48) 



The high-frequency spectrum is thus purely sinusoidal and has no modes or divergences. However 
the location of the minima and maxima, their separation, and their visibility are dependent on 
both the type and location of the source. 



3.3. The Intermediate-Frequency Case 



Formally we can consider a third case in which the upper turning point lies well above the 
source. In practice this case probably does not accurately describe the properties of any actual 
mode occurring in the Sun, so I include it only for completeness and note only the result for the 
amplitude spectrum: 



F 



1 



2 cos ( 



ao sin 



7T 

+ 4 



ai cos 



7T 

+ 4 



where <fi is given by equation (29) and 



Kdz. 



(49) 



(50) 



The condition for the eigenmodes is the same as in the low-frequency case [eq. [3l| . However the 
absence of the small factor exp — 2A means that the spectral peaks are much more symmetrical in 
this case. 
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4. Application to a Simple Model 

Asymptotic descriptions of the properties of solar oscillations have generally proved useful in 
helioseismology as guides to the relationship between properties of the spectrum and properties of 
the model. Only in limited cases can the actual asymptotic expressions themselves be used for 
direct numerical comparison between observation and theory. In the present case we have found 
that the asymptotic analysis points to the phase-function difference between the modes and the 
troughs in the solar spectrum as being a useful diagnostic for the source of the oscillations. I now 
apply this to a simple model. 

I consider a model consisting of an adiabatic polytropic layer characterized by a polytropic 
index m, truncated at some depth z s corresponding to the surface of the Sun, and overlain by an 
isothermal atmosphere characterized by a density scale-height, -£/_. A parameter (3 is defined as 
the ratio of the temperature immediately above the surface to that immediately below. Thus for 
z > z s we have 

4 = il H + = - u 2 c+ = ^-(l + -) and N 2 + = (51) 



m m 4z \ m J 



while for z < z s we have 

2 o9z s tt n z s 2 lg(m + l) 2 2 g 

c_ = (3 — H- = P—-r u c _ = -— and N_ = — . (52) 

m m + 1 p 4mz s pz s 

Figure [5] shows asymptotic power spectra for an m = 3 polytrope in which a pure monopole 
source is located 90km below the photosphere. The peaks show the same sign of asymmetry 
as seen in the solar doppler observations. The asymmetry decreases with increasing horizontal 
wavenumber and increasing frequency. Also shown is the effective difference between the phase 
function describing the peaks and the phase function for the troughs. Evidently the difference 
is close to being a function of frequency alone, even towards the higher frequencies where the 
approximate expression for the difference [eq. |3^] is invalid. 

Figure [6| shows the spectrum for the same case as in figure || however I now show the result 
for a single wavenumber but extended over a large frequency range so as to include all three 
asymptotic regimes. The crossovers between the three asymptotic regimes are quite smooth despite 
the fact that the asymptotic analysis formally breaks down at those points. The asymmetry is 
again strongest at low frequency and decreases as the frequency is increased. However there is no 
actual reversal of the sign of asymmetry for this particular source depth. In the high-frequency 
range the troughs are still present but there are no eigenmode peaks. When the monopole source 
is replaced by a dipole the locations of the eigenmodes are, of course, unaffected but the sign of 
the line asymmetries is reversed. For the high-frequency modes there is a phase shift so that the 
locations of maxima and minima in the spectrum are interchanged. 
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4.1. Comparison of Asymptotic and Numerical Results 

For comparison between the asymptotic and numerical results, I solve equation (||) with 
right-hand-side given by equation (^2) using a shooting method. Specifically one starts from 
a location well above both the upper turning point and the source with the initial condition 
\P' — \K\ty = and integrates down to the source at zq. Call the values of the wave function and 
its derivative there \P_ and ^>'_. Similarly one can integrate from below the lower turning point up 
through the cavity to zq, starting with the boundary condition + \K\^f = 0, to obtain ^_|_ and 

immediately below the source. Then the jump conditions [eqs. |1| 16 imply that the actual 
solutions above and below the source are given by j4_^_ and A + ^ + where the coefficients A_ 
and A + must satisfy 

A+V' + - A^'_ = aol-fifC^)^ (53) 

and 

A+V+-A-ii!- = a 1 \K(zo)\-* (54) 



giving 

The condition for an eigenmode is therefore that the denominator vanishes, while the condition 
for a zero in the power is that the numerator vanishes. Without loss of generality we may assume 
that fy+ and are real. If we consider only cases where ao and a\ are in phase, or cases in which 
either ao or a\ vanishes, we see that zeros of A- can occur and that their location is independent 
of Thus the location of the nulls in the spectrum is independent of the atmospheric structure 
above the source. The explanation for this result is quite simple - since the nulls correspond to 
frequencies where the amplitude in the superficial layers is identically zero their location cannot 
depend on the structure of those layers. In the more general case in which ao and ai are complex 
the minima will no longer be true zeros and their location will no longer be independent of 
Instead they will be given by the condition d|A_| 2 /dw = which in general will depend also 
on Nevertheless, the minima will still be locations in which the amplitude in the region 
above the source is small and their location will therefore still be less sensitive to the structure 
of the superficial layers than the locations of the eigenmodes. This result is illustrated in figure 
which shows a single peak and trough in the case of a monopole and dipole source acting in 
phase-quadrature and for two different values of the photospheric temperature-jump parameter. 
Evidently the trough remains almost fixed while the eigenfrequency varies. 

If we return to the asymptotic expressions for the eigenfrequencies and nulls [eqs. |3l] and 



32 1 we see that these expressions imply that the mode frequencies depend only on the structure 
below the upper turning point while the nulls depend only on the structure below the source. 
However, while the former statement is only an approximation, because the eigenfunctions have 
non-zero amplitude outside the turning points, the latter statement is exact for those cases where 
the troughs are true nulls. 
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Figure [8| shows a comparison of numerical and asymptotic estimates of the power spectrum. 
While the spectra are qualitatively the same, quantitatively the separation between the peaks and 
troughs for the numerical calculation is much smaller than that predicted by the asymptotics. 
This suggests that while the asymptotics may be a useful guide as to the correct interpretation of 
the line asymmetries, it would be wrong to use an asymptotic expression such as equation (35) in 
a quantitative comparison with the data. 



5. Discussion 

Typically, helioseismology has concentrated on analysis of the eigenmode frequencies of the 
Sun and relatively little attention has been paid to the rest of the spectrum. With the advent of 
the SOI/MDI instrument on SOHO and the GONG network, noise levels in measured spectra are 
now sufficiently low that it is not only possible to make useful measurements between the mode 
peaks but necessary to do so since it is increasingly clear that the line asymmetry can produce 
systematic errors in the determination of p-mode frequencies if it is not accounted for properly. 

A puzzle in the data has been the observation originally due to Duvall et al (1993), and now 
apparently confirmed from SOI-MDI (Nigam et al. 1998), that the sign of the measured asymmetry 
is oppositely directed in doppler and intensity measurements. The doppler measurements imply 
a dominant monopole source (a negative asymmetry in the Abrams & Kumar terminology) while 
the intensity measurements would suggest a dipole source. Rast & Bogdan (1997) show that, 
regardless of the form of the source function, one would expect to see the same sign of asymmetry 
with observations made in different variables, since the locations of the peaks and troughs are the 
same for all variables. Nigam et al. have shown that the phenomenon can be explained if the 
measured signal in either doppler or intensity is contaminated by a noise source which is correlated 
with the source, for example if the measured signal contains a contribution coming directly from 
the source itself. The question is then raised as to which of the doppler measurement and the 
intensity measurement actually constitutes a better measurement of the true acoustic spectrum. 
One might argue that the doppler signal is probably less contaminated since it has a lower overall 
background level, but it is certainly possible that both the intensity and doppler measurements 
contain some leakage from the source, with an amplitude and phase which are not yet precisely 
known. Evidently the possibility that the location of the spectral minima is shifted due to such 
leakage must complicate the interpretation of the line-asymmetry in terms of source location and 
type, and further work will be required to determine whether the combined doppler and intensity 
measurements can be used to determine the true acoustic spectrum in a model-independent 
fashion. 

As originally pointed out by Rast & Bogdan, the asymmetry is best approached by 
consideration of the minima in the spectrum, whose location is determined by the depth and type 
of the source of the oscillations. In this work I have determined an asymptotic expression for 
the frequency spectrum, from which it is possible to determine the locations of the minima. The 
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sign of the line asymmetry is determined by whether the nearest minimum lies above or below 
the eigenfrequency. The former case corresponds (at low frequency) to a pure monopole source 
and the latter to a pure dipole. The amount of asymmetry depends of the location of the source. 
The higher the source is placed above the upper turning point of the oscillation the closer the 
trough will lie to the eigenmode and therefore the greater the asymmetry. Thus the extent of the 
asymmetry should be a useful diagnostic of the properties of the turbulent superadiabatic layer in 
which the source of the solar p-modes is believed to lie. 

Numerical results confirm the qualitative properties of the asymptotic results. Also for the 
case of a pure monopole or dipole, or where the two components of the source are in phase, the 
exact locations of the nulls (including those above the acoustic cutoff frequency) can be determined 
by solving numerically the oscillation equations subject to the boundary condition 

^-^V|* = (56) 

at z = zq. This should be relatively straightforward to implement in most existing oscillation 
codes. One can thus determine numerically the phase-function difference at — a m which is 
predicted by the asymptotic analysis to be a key diagnostic of the structure and location of the 
source. 

Since the nulls are independent of the superficial structure the asymmetry could, in principle, 
provide a tool for differentiating between superficial and more-deeply seated changes in the Sun. 
In particular, solar-cycle variations occurring above the source location, in the lower atmosphere 
for example, would change the mode frequencies but not the trough frequencies and could thus be 
visible as changes in the asymmetry, while changes below the source would change the frequency 
of both the mode and the trough. In practice such changes are likely to be extremely difficult to 
observe, not least because the asymmetry is most evident for lower frequency modes for which the 
solar-cycle dependent changes are small. 

We have thus reached the stage where it is possible to envisage a program for the re-analysis 
of helioseismic data to make full use of the information contained in the p-mode line profiles. 
In particular it would be highly desirable to re-analyse the spectral data using parameterised 
line-profiles such as equation (|42|) or equation fl44| ) in order to obtain improved estimates of 
the eigenfrequencies simultaneously with estimates for the spacing between the maxima and 
minima. From these we could obtain an observational at — a m curve which could be calibrated 
against theory using nulls calculated from an eigenvalue calculation with boundary condition ([56]) . 
Ultimately one would hope that the results of such a calibration could be brought into line with 
predictions made from numerical simulations of the solar surface layers (e.g. Stein & Nordlund 
1989, Stein & Nordlund 1991, Nordlund, Stein & Brandenburg 1996) to provide a solidly grounded 
understanding of this highly complex region of the Sun. 
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Fig. 1. — A schematic of the illustrative model used in 



Fig. 2. — Power (solid line) and energy (dotted line) spectra for the illustrative model for three 
different sets of parameter values. 



Fig. 3. — A schematic illustration of the geometry for the low-frequency regime analysed in §3.1 



Fig. 4. — Preliminary results of a least-squares fit of equation ( f42|) to m-averaged power spectra 
obtained from one month of GONG data. Each panel is an average over all p-mode peaks in the 
range < / < 10 in the frequency bin indicated, with the mode peaks aligned using the most 
recent GONG tables of frequencies and splitting coefficients. The dashed line shows the fit to 
the asymmetric profile (^) and the dotted line is a fit to a Lorentzian profile with a constant 
background. 



Fig. 5. — Asymptotic power spectra resulting from driving an m = 3 polytrope with a monopole 
source located 90km below the photosphere. Results are shown for three different values of the 
mode degree I, defined by I = kR & . The lower panel shows the difference in phase function between 
the trough and the mode for each trough-mode pair appearing, using the symbols (o) for I = 100, 
(+) for I = 200 and (*) for I = 300. 



Fig. 6. — Asymptotic power spectra over a broad frequency range for the same polytropic index 
and source depth as used in figure || for I = 200. The solid line is for a pure monopole source and 
the dotted line is for a pure dipole source. The crossovers between the three asymptotic regimes 
are marked by vertical dashed lines. 



Fig. 7. — Numerical power spectra for a single mode with I = 300 for the two cases (3 = 1.5 
(dotted line) and f3 = 10 -3 (solid line) for a source characterised by ao/a% = — lOi (and with other 
parameter values as in figure |||) . 



Fig. 8. — Comparison of numerical (dotted line) and asymptotic (solid line) low-frequency spectra 
for I = 200. The lower panel shows the frequency difference between a trough and its neighbouring 
peak for each pair in the calculation. 



